SCATTERING FOR SOLUTIONS OF NLS IN THE EXTERIOR OF A 2D 

STAR-SHAPED OBSTACLE 



FABRICE PLANCHON AND LUIS VEGA 

Abstract. We prove that solutions to non-linear Schrodinger equations in two dimensions 
and in the exterior of a bounded and smooth star-shaped obstacle scatter in the energy 
space. The non-linear potential is defocusing and grows at least as the quintic power. 



1. Introduction and main results. 
In this paper we are interested in the initial value problem 

{idtu + Au = e|n|P"^n n G \ S = , e G {0, 1}, p > 5, 
u\g^ = 
u(x,0) = uo{x)eH^{n) 

and more particularly in the proof of scattering properties of the solutions. With respect to 
E we will assume that it is a star-shaped domain contained in a compact set K, and that 
its boundary dH = dQ is a smooth curve in M^. 

As it is well-known dimensions one and two are the most delicate ones regarding scatter- 
ing questions even if no obstacle is considered and one works in the full space M'^. The 
main obstruction comes from the sign of the bilaplacian term that appears in the classical 
Morawetz-Lin-Strauss inequality [9]. This sign turns out to be the wrong one for d = 1,2. 
In dimensions three and higher and in M*^ the question was settled in the energy space by 
J. Ginibre and G. Velo in |6j . The obstruction in the dimension was finally removed, still 
in full space, by K. Nakanishi, [10], but his technique does not seem directly applicable in 
the exterior of a 2D domain; that is our interest here. One should remark that the ID 
case (with both Dirichlet or Neuman boundary conditions) follows directly from the M case, 
as symmetry with respect to x = is preserved by the nonlinear flow. Hence for exterior 
domains the 2D case may be seen as the most difficult one. 

Another fundamental contribution regarding global existence and scattering was the in- 
troduction by Colliander, Keel, Staffilani, Takaoka, and Tao in |3] and |5] of interactive 
Morawetz inequalities. However a similar obstruction in low dimensions appears as well. 
Again the difficulty comes from the sign of the bilaplacian term that appears in the use of 
the bilinear multiplier. It turns out that this obstruction can be easily bypassed (see [2] and 
and as a consequence the scattering of NLS equations for subcritical non-linearities 
is by now fairly well understood (see for example [7j for a good survey, and also [3] for earlier 
results in ID using a quadrilinear interactive Morawetz inequality). 
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In [H] the bilinear multiplier technique is also used in the setting of exterior domains with 
Dirichlet boundary conditions. In order to explain our results we have to introduce some 
notation. 

Let 

n>l , p>l , eG {-1,0,1} , Qi,Q2CW 
with smooth boundaries dQi,dQ2, and Ui,U2 solutions of 

(2) idtUj + Auj = e \ujf~^ Uj with WjlgQ. =0 j = 1,2. 

Denote for j = 1, 2 

(3) M{uj)= I \ujfdx; E{uj) = ]- j \Vuf + f \ujf^^ dx, 

Jqj 2 Jq^ p + 1 Jq. 

the mass and energy that are conserved quantities. As in [TT] we shall use the Radon 
transform 



(4) R{f)is,u)= / /d/i,,^ 

J {xuj=s}nfi 

with Q either Qi or Q2- In other words we extend / as zero outside of Q and compute its 
Radon transform. 

We set 

(5) Ip{t) = Ip[ui{t),U2(t)) = p{x - y)\ui{x)\'^\u2{y)\'^dxdy. 

Then a simple modification of the proof of Theorem 2.5 in |11] gives us0 

Theorem 1. Let uj G M", n > 1, with \uj\ = 1 and Pu]{z) = \z ■ uj\, Uj solution of ([2]) with 
j = 1,2. Then if x = x-^ + su and x^ ■ u = we have 

(6) ^^Ip^= j\ds{n{u,U2))fds 

7^(|nl|')7^(|n2|^+')rfs+ / 7^(|u2|')7^(|nl|P+l)ds 



p — 1 



s J x-bj=s J y-u)=s 



1 [x^ + suj) dsU2 [y^ + sbj) 

U2 {y^ + su) dgUi (x"*" + su) 1^ dx^dy^ds 
U2\'^{y)dnpLo{x - y) \dnUi\^ {x)dS^dy 



\ui\^{x)dnpui{x - y) \dnU2\^ {y)dxdSy 



where dn is the outgoing normal vector field on dQ. 



^There are some misprints in the second and fourth term of the RHS of (2.19) in [TT]. A factor 2 is missing 
in both of them. 
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A simple look at the statement of the above theorem tells us that it is useless unless the 
boundary terms appearing in ([6]) are under control. In [11] this is done using the so- 
called local smoothing property that follows from a variation of Morawetz's multiplier. As 
a consequence a sufficient geometric condition on the boundary naturally appears: if u is a 
solution to ([2]), and n{x) denotes the outgoing normal vector to dQ at x, local smoothing 
for the nonlinear solution u is obtained provided 

[ [ {x ■ n{x))\dnu\^ dS^dt > . 
Jo Jdn 

In particular, it is sufficient to assume that Q is the complement of a star-shaped domain 
and this is what we are going to assume in this paper. Unfortunately the use of Morawetz's 
multiplier introduces again the restriction on the dimension and therefore the result in [11] 
(Proposition 2.7) is given in dimension three and higher. In this paper we remove such a 
restriction. 

The main new idea is to use Theorem [1] with ^2 = f2i = Ui = u and U2{x) = u[—x) and 

u a solution of ([T]). Equivalently we may consider p = p{x + y) instead of p{x — y) in the 
definition of Ip in We then sum up the new ( with p^) and old (with p^) estimates: it 
follows that some cancellation occurs between boundary terms of and integrating over 
u G one obtains that the boundary term is bounded above by 

/ \dnu{x)\'^\u{y)\'^{y)~^dS^dy , 

JdnxQ 

1 /2 

where the gain in the above inequality comes from the weight (y) = (l + 

It turns out this gain is sufficient. The reason is that we can use the tensor product technique 
developed in [^ and |3]. In order to do that we construct v{x,y) = u{x)u{y), solution of 
NLS in f2 X f2, and use the local smoothing inequality obtained from Morawetz's multipliers 
in dimension n = 4, see Proposition [1] below. That fits perfectly well with our purposes. 

Remark 1. In the case, all monotonicity formulae stemming from virial identities are 
related to the conservation of the momentum, a key property which is lost on domains. 
The boundary terms then appear as a direct consequence. The star-shaped case provides a 
class of obstacles for which this loss is controlable by an integration by parts argument (the 
usual virial in dimension n > 3 or the tensorialized one in the present paper for n = 2). 
For the linear equation, microlocal techniques allow to control the momentum through local 
smoothing estimates - as they control boundary terms, see [11] when the exterior domain 
is non trapping (light rays escape to spatial infinity). But extending such a control to the 
nonlinear setting appears a challenging task. 

We denote by = Hq{Q) the energy space, which is also the domain of the square root of 
—Aq, where Aq is the Dirichlet Laplacian on Q. We are now able to state our main theorem. 

Theorem 2. Let Q beM.'^\T,, where is a star-shaped and bounded domain, anduo G Hq[Q). 
Then, there exits a unique solution of 

\u\p~^u xen , teR p>5, 
0, 

uo{x), 



idfU + Au 



u 



an 



u{x, 0) 
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such that 



ueC{R:H',{n))nLr'L:^. 



Moreover, there exist unique u+,u^ G Hq{Q) such that 

(7) -e^*^"n±|| = o(l) t ^ ±00. 

Here e**^" is the solution of the linear problem (i.e. ([1]) with e = 0). 

The reader can guess that the key property of the solution obtained in the above theorem is 
that it satisfies the global in time Strichartz estimate L^~^L^ . The key step in that direction 
stems from Theorem [T] and the local smoothing estimate that we mentioned before: one gets 
that is in of both variables space and time. This is also true for linear solutions. 

This turns out to be the last key ingredient for the proof of Theorem |5J We have the following 
result. 

Theorem 3. Let f2 = \ S, where S is star-shaped and hounded. Define e**^"uo as the 
solution of 

' idtU + Am = X en , t eR 



Then 

and, for < s < 3/4, 
As a consequence 



u{x,0) = uo{x), 
\\D^u\\,,^s<C\\uo\\^s, 



-1/4 



Here = (— A)'^/^ and Hq{VI) = Hq denotes the domain of the operator (—AqY^'^ given by 
the spectral theorem applied to the Dirichlet Laplacian on Q. 

Remark 2. Notice that under the smoothness conditions on Q, (— A^)"^^^ is hounded in 
1 < p < 00. See for example |8] and references there in for the proof of this fact. 

Remark 3. One may extend Theorem\^ to nonlinearities with p < 5, hy using the known 
range of Strichartz estimates outside non trapping obstacles from [I].- for e > 0, 



\u\ 



<C(e)||no|| 



This allows to replace p = A in TheoremlM by p = 3/{l — e), which will however not fill 
the expected range 3 < p < 5 in Theorem\^ by analogy with the case, but rather provide 
p > 4. We elected to keep the argument mostly self-contained, as the linear estimates from 
Theorem\^ are enough for our purposes in the range p > 5. 

In the next section we prove Theorems [2] and [3] and give the necessary propositions to be 
able to use Theorem [H This one is proved as Theorem 2.5 in and therefore we omit the 
details. 



scattering for nls outside a 2d star-shaped obstacle 

2. Proofs. 

We will first prove Theorem [3l We need some lemmas. 
Lemma 1. For u G define Puj{x) = \x ■ uj\ and 

^ui = {{x, y) e dil X : \y ■ u\ < \x ■ u\} . 

Then if y E Q 

{{x,y) e on X n : dnpuj{x -y) + dnpu^{x + y) 7^ 0} 
C {sig {x - y) ■ u + sig{x + y) ■ uj ^ 0} = A^, 
and 

C 



(9) 



(10) 



sup 

xedn 



< 



(y)' 



Proof. We can assume without loss of generality that u = (1,0). Then if x = {xi,X2), 
y = {yi, 2/2) we have that 

dnPuj{x - y) + dnPuix + y) = sig{x - y) ■ u + sig{x + y) ■ u 

= sig (xi - yi)) + sig (xi + yi)) 

and hence by inspection we get ([9]) 

As for ( JTOj) we have that if ||/ ■ u;| < \x ■ u)\ then \y\ \cos{y,uj)\ < diam (dQ). Thus, for \y\ 
large enough 



9y — 9^ ± — 



< M 



with ill = (cos^y, sin6'y) and u = {cos6^,sm6^). The lemma easily follows. 



□ 



Proposition 1. Let Q 6e \ S, with S a bounded smooth star-shaped obstacle. Assume u 
is a solution of 

idtU + Au = e|M|^~^M, p > 1, 



u 



an 



0, 



with e G {0, 1}. Then 
ill) 



I I \dr.u{x)\'\u{y)\'dS.jLdt<CM^''E'l' 
jRJdnxn {y) 

with M and E as in 

Proof. Define v{x,y) = u{x)u{y) solution of the problem 

idtV + Av = e{\u\P-'^{x) + uP-'^{y))v, (x, y) G x fi, 



(12) 



\d(nxn) 



0, 



v{x,y,0) = uo{x)uo{y). 



Then, consider for h{x, y) = ^/\x\^~+~\y^ 



\v\'^{x, y, t)h{x, y)dxdy 



nxn 
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and compute -j-^Mh(t). The details can be found for example in (p. 278), up to easy 

modifications to handle the nonlinear term in ( fT2l) . We just have to bother about dnh with 
n the normal to d{Q x Q), namely 

n = {rix, 0) if X G dQ , y ^ Q, 
n = {0,ny) if X E dQ , y E Q. 

Hence 

dnh{x, y) 



if (x, y) G dVL x f2 and similarly in the other case. We recall that the kernel of a star-shaped 
domain is the set of points with respect to which the domain is star-shaped. If we assume 
that the kernel of our obstacle contains a disk (note that the kernel is always convex) then we 
may exclude the situation where n^. - x = 0, by averaging over this disk if necessary. Abusing 
notation by forgetting about this possible average over the base point, we then get that 

C 

\dnh{x,y)\ > — 

\y\ 

with C > 0. The argument follows as in pg 278 of [H] because 

\dnv\^ = \dnu{x)f \u{y)f 
if (x, y) G dQ X Q and a similar expression if (x, y) G x dQ. □ 

If the kernel of the obstacle is a single point (or even a segment), we need to address 
the situation where there exist points x G dQ such that n^. ■ x = 0. We may repeat the 
computation of [H] with weight p(x, y) = a/1 + |xp + but the desired control on the 
boundary does not follow by a subsequent computation with a weight like the distance to 
the boundary, because in our setting the 4D obstacle is not compact and we need to account 
for the trace of dnh on the boundary. Let us summarize what the computation with weight 
p provides, following 

Lemma 2. In the conditions of PropositionUl 

(13) // ^ ^ IVH; + I'f ^ (I'-l'-M + \ur'(y))H' < _ 



jRJnxnp{x,y) P^{x,y) pix,y) 
where ^ denotes the (4D) angular gradient. 

We now proceed with the remaining case in the proof of Proposition [1] : control of the 
boundary term in the situation where there exist points x G dQ such that x ■ n(x) = 0. 
While Lemma [2] is not enough to conclude directly, it will be a key ingredient in what 
follows. Recall we aim at controling the lefthand side of ( ITT]) . 

Let R be such that dn C B{0,R-1). For the part of the 4D boundary {dnn B{0,R)y we 
may proceed as in using the control of Vf in f|T3|) . because the weight is irrelevant. 
The remaining part of the 4D boundary is a union of cylinders dQ x (M^ \ B{0,R)) and 
(M^ \ -8(0, R)) X dQ. They will be treated similarly and we restrict our attention to the first 
one (boundary in x and exterior of a large ball in y). We may further divide the region in 
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a finite number of angular sectors: witliout loss of generality we now assume to be in such 
a sector 5* = {\y2\ < |yi|/10}- Notice how being outside a ball forces \yi\ > R/2. 

From the compactness of dQ and its smoothness as a curve in M^, there exists a vector field 
Z such that Z is defined in Q, its restriction on dQ is the normal derivative, and it vanishes 
outside B{0,R). In a coordinate system, 

Z = a{x)d^^ + b(x)d^^, 

where a{x), b{x) are smooth functions vanishing for \x\ > R. We define two 4D vector fields 
on our 4D region of interest, {Q \B{0,R)) x S: 

Yk = Vidxk - Xkdy^ , with k = 1,2. 

Notice that Yi and Y2 are, up to a factor (a/I^P + lyP)"""", angular derivatives (in 4D). We 
may now define 

Y = a{x,y)Y, + f3{x,y)Y2 

where the support of a and /3 is in {Q\ B{0, R)) x 5* and their respective value is adjusted 
such that (Yv)^xedn ~ h~^{x)Zv, as dy^^v = on dQ x Q: set a{x,y) = 4>{y /\y\)(i{x)/yf 
and (]{x,y) = (f){y /\y\)b{x)/yf where (j){y/\y\) is the identity in the angular sector S, and 
vanishes outside 2S. One should think of the following computation as an analog of the one 
in [H] p. 278/279, but where we think of the region of 4D space we are integrating over as 
being foliated by spheres |xp + |yp = C rather than planes \y\ = C. We use the following 
momentum 



J = 2Im / v{Yv) dxdy. 
Jnxn 



For the sake of notational convenience, we set X = [x, y) and Y = A[X) ■ Vx = A ■ V and 
N{u) = \u\^~^{x) + \u\'^~^{y). The following computation is standard, does not depend on 
the specific form of A and similar to the one in [11] where A = V/i; 



dt ^2Im j vVv ■A^=-2lmJ dtv{2Vv ■ A + vV ■ A) 



-4Re j /\vVv-A + 2 j \Vv\'V -A- / ■ A 

2 I NV{\v\^) ■ A+\v\^NV -A. 



We have then (with An the normal component of A on the boundary and da the surface 
measure on it, with orientation toward the domain) 

2Re j AvVv ■ A = -2 J \dnv\'^Anda - 2Re J Vv ■ V(A ■ Vv) 

= - 2 y \dnv\'^Anda -2Re j A- V(Vw ■ Vv) 

+ 2Re y ■ V, V]v ■ Vv 
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and by substitution and integration by parts in the middle term and the nonhnear term 
above. 



dt (^2lm J vVv -Aj = 2 j \dnv\'^Anda - 4Re ^ " V, V]i; ■ Wv 



- I \v\'AV-A + 2^ I \v\\\ur\x)V,-A, + \ur\y)VyAy), 



V + 1 

where A^ (resp. Ay) stands for the projection over the x coordinates (resp. the y coordinates) 
of the vector field A. Given our choice of F = A ■ V, the boundary term is what we seek, up 
to an harmless factor h{x,y)/yi: 

f \dnv\^Anda= [YvVvda= [ tMMl\d^y\^ dS^dy . 
J J Janxn Vi 

The term N{u)\v\'^ will be under control by Lemma [2l For the |fp term, we may use 



2 



Poincare's inequality on the spheres |a;p + = C to control it as well by J p ^(x, y)\^ 
(notice how when deriving in x we do not gain decay in the coefficients of Y). 

We are left with the term carrying [y,V]: for any scalar function 7, we have [7yfc,9] = 
7[yfc, d] — {d'-f)Yk = —{d'~f)Yi: ; hence the vector field [Y, V] is spanned by angular derivatives. 
Moreover, \d'~f\ < I7I (as the worst situation is when the derivative would hit a(x) or b{x), 
which does not gain decay in y) and 7 ~ yi/p{x,y), so that 



\[Y,V]v-Vv\ 



< 



2 



pix,y) 

which is controlled by Lemma [2] after time integration. This concludes the proof of Propo- 
sition [U □ 

Proposition 2. In the conditions of Proposition^ 



(14) 



/ / \D'/^\u\''\'dxdy<CM'/^E'/\ 
Jr Jnxn 



Proof. We use Theorem [T] with Ui = u{x) and U2 = u{—x), Q' = —Q. Then if y & Q' 

dnPu{x -y) = dnPu{x + y). 

We use again Theorem [T] with u = Ui = U2 and f2i = ^2 = ^- Then we sum the two 
identities and integrate in time once. After discarding some positive terms we get 



T 







|9s7^|np(s, oj)^ dsdu 



2 / \dnu{x)\^ \u{y)\^ {dnPu{x -y) + dnPu{x + y)) dS^dy 



anxn 

2 / \dnu{y)\^ \u{x)f {dnPuj{x - y) + dnPuj{x + y)) dxdSy 
Jnxdn 



where the right hand side follows from the trivial bounds of /^^ 
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We conclude the argument integrating on a; G and using the well known Plancherel's 
Theorem for the Radon transform as in (Prop 2.2) to obtain 



1^1/2(1 



U\ 



|2\ 



r 2 

x,lO,T] 



For the boundary term we use Lemma [T] and then Proposition [H Notice that both lemmas 
are symmetric in the (x, y) variables. □ 

Remark 4. The astute reader will have noticed by now that the entire computation may be 
performed directly, with the (integrated over u) weight p{x,y) = \x — y\ + \x + y\. We feel 
the directional weight to be more natural, with the key LemmaUl having a simple geometrical 
proof. In fact, we were led to the present result from considering simple obstacles with 
symmetries (e.g. the disk). 

Proposition 3. In the conditions of LemmaU\ and e = 0, 

(15) WD'^uW^^s < C\\uo\\^s,+i/4 

Proof. Take uq a frequency localized function (according to the Dirichlet Laplacian of Q). 
Then ffTSj) is a consequence of (fl^ and Sobolev embedding. The general case follows by 
the square function estimates on domains and the boundedness in L^, 1 < p < oo, of the 
corresponding Riesz transforms, see for example j8], Theorem 1.2 and Remark (1.10). 

□ 

Notice that Proposition |2] and Proposition [3] imply Theorem [31 We already have all the 
necessary ingredients for the proof of Theorem [2l 

Proof of Theorem\^ For simplicity we will consider the particular case p = 5. The necessary 
changes needed for the rest of the powers will be outlined at the end of the proof. 

We first notice that by Gagliardo-Nirenberg inequality 
Hence 

(17) l|e'*^"«o|L4^c, < C\\u4h^^, s > 1/2. 

This inequality together with the trivial one 

u{t) = e^*^«no + I [ e^(*~")^"|n|^n(r)dr. 



IVImI'^mILo < c||m||^oo II Vm||l2, 
allows us to do a fixed point argument in L^L'^ fl L'^Hq, to solve locally in time 



^0 

On the other hand given any interval / of time where the solution exists we know from 
Proposition [2] that 11^11^42^8 is finite and bounded by E^^^M'^^^. Hence given any e > there 

N 

is a finite number of disjoint intervals Ii, . . . ,1^ such that [J Ij = I with = A^(e) and 
IMlj L°° = ^' J < ^(e)' IMlj L°° ^ J = N{e). 
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H6I1C6 

Assume Ij = [tj,tj+i) so that 

\Ht,)\\H^ < E'/\ 
and call ^ 



Then a{t) is a continuous function with critj) = 0. Moreover for tj < t < tj^i 

ll«''*''L?,,.|£C'^"'^(l + ll-lllf,,.r)- 

Hence 

< Ce2/3E2/3(l + a(t)4/3)^ 

and taking eE small enough, (i.e. e = {CE)~^) we conclude that a{tj^i) remains bounded 
by a universal constant independent of /. Therefore the solution is global and 

||M||r4.oo <CN < CE^/'^M^/\ 

Finally defining 

/■oo 

u^ = Uo+ / e'^^''\u\%{T)dT 
Jo 

we obtain ([7]). 

For the general case p > 5 we need a bound for ||n||^p-i^^. Call q = p — 1. Then for 6 = 4/q 
we have 

"+00 



^x 

Hence 



/" / \ 1 

'n\\i,a^<c (iklllsikii^/) dt<cM-'i^Ei-^. 

'J 



Finally it is enough to use instead of ( !T6|) the following inequality 

The rest of the proof works similarly. □ 

Remark 5. Notice how, in the previous argument, one is combining a Strichartz estimate of 
type LIL^ with the norm to control the source term \u\^~'^u in L\{Hl), and this in turn 
yields r = p—1. Given that all available Strichartz estimates on domains are suboptimal with 
respect to regularity - there is a loss in the Sobolev embedding scale when compared to - we 
cannot hope to use e.g. dual Strichartz pairs to allow for higher time integrability . As such, 
dealing with a given p requires the availability of Strichartz estimates with r = p — 1. In light 
of ongoing developpement on linear dispersion on domains, we may hope (being optimistic) 
for r > 12/5, which is still not the range, and new insights on the linear theory will be 
needed to match it. 
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